Dynamics of rough interfaces in Chemical Vapor Deposition: experiments and model 

for silica films 
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We study the surface dynamics of silica films grown by low pressure chemical vapor deposition. 
Atomic force microscopy measurements show that the surface reaches a scale invariant stationary 
state compatible with the Kardar-Parisi-Zhang (KPZ) equation in three dimensions. At intermediate 
times the surface undergoes an unstable transient due to shadowing effects. By varying growth 
conditions and using spectroscopic techniques, we determine the physical origin of KPZ scaling to 
be a low value of the surface sticking probability, related to the surface concentration of reactive 
groups. We propose a stochastic equation that describes the qualitative behavior of our experimental 
system. 
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The dynamics of growing interfaces has become of in- 
creasing interest in order to understand the physical pro- 
cesses that determine film quality Q. In the absence of 
morphological instabilities, many surfaces evolving out of 
equilibrium display time and space scale invariance, i.e. 
such surfaces are rough. A successful framework for the 
study of rough interface dynamics has been the formu- 
lation of stochastic differential equations for the surface 
height h(r, t), where r denotes a site of a two-dimensional 
substrate, and t is time. A prominent example is the 
Kardar-Parisi-Zhang (KPZ) equation 



dth = vV 2 h + (A/2)(V/i) 2 + ri(r, t), 



(1) 



expected to describe the large length scale dynamics of 
any rough surface growing in the absence of specific con- 
servation laws. In Eq. (|l|), v and A are constants. The 
first term on the rhs is sometimes called a surface ten- 
sion term, as it tends to smooth the interface through 
evaporation-condensation processes Q. The nonlinear 
term accounts for growth along the local normal direc- 
tion (lateral growth), and T)(r, t) is a Gaussian white noise 
of constant strength 2D, accounting for microscopic fluc- 
tuations, e.g., in a deposition beam. 

One of the techniques for thin film production which 
might be expected to lead to surfaces described by (Q) is 
Chemical Vapor Deposition (CVD). CVD is very widely 
used for technological applications and industrial devices 
||. It is also very interesting from the fundamental 
point of view due to its conceptual similarities with other 
growth techniques such as electrodeposition. Thus, a 
considerable effort has been made to model CVD growth 
||-|| using partial differential equations. These models 
predict the development of unstable morphologies within 
the time ranges studied, contradicting the above expec- 
tation of KPZ scale invariant behavior. Let us note 



that, to date, very few experimental systems are known 
whose scaling is described by KPZ in three dimensions 
[Q. Moreover, on the fundamental level there are no ex- 
periments addressing the long time behavior of the sur- 
face morphologies produced by CVD, and the identifica- 
tion of the determining physico-chemical mechanisms. 

In this Letter we study the interface dynamics of low 
pressure Si02 CVD films by Atomic Force Microscopy 
(AFM), Raman and infrared spectroscopies, and numer- 
ical simulations. Si02 films were studied because of their 
amorphous nature, to avoid Schwoebel barrier effects on 
the surface dynamics and to prevent formation of 
facets that can also alter the scaling behavior Our 
experimental system displays a transient unstable behav- 
ior related to the gas phase transport character of CVD 
growth Q-^J, but a stationary state is achieved which, 
rather, is compatible with KPZ scaling. The surface dy- 
namics is well described, over a wide spatial and temporal 
range (up to 2 days of deposition time), by a stochastic 
differential equation in which surface diffusion, shadow- 
ing effects, and lateral growth are allowed for. 

Amorphous Si02 films were grown at a temperature T 
— 723 K on Si (100) substrates in a hot wall, horizontal, 
low pressure tubular CVD reactor. The precursor gases 
were silane (diluted at 2% in nitrogen, 99.999% purity) 
and oxygen (99.9992% purity), with an oxygen/silane ra- 
tio equal to 20 and a total gas flow rate of 50 seem. The 
chamber pressure was 1.4 Torr. The film thickness in- 
creases linearly with deposition time at a constant growth 
rate of 20 nm/min. Films were deposited in the range 5 
min < t < 2 days. The surface morphology was char- 
acterized by AFM (Nanoscope III from Digital Instru- 
ments, CA) operating in tapping mode at ambient con- 
ditions up to a scale of 50 /im using silicon cantilevers. 

AFM imaging of our silica films shows (Fig. 1) a de- 
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posit formed initially by small rounded grains 30-60 nm 
in size. A similar morphology has been reported on amor- 
phous silicon films deposited by thermal evaporation |j] . 
As deposition proceeds, structures resembling mountains 
and valleys appear at larger length scales increasing in 
size, until a stationary regime is attained in the 15-30 
hours range. 

Power Spectral Density (PSD) (or surface power spec- 
trum) plots of the surface morphology at different times 
are displayed in Fig. 2. For t < 35 min and small enough 
wave vector k, the PSD follows that of the initial sub- 
strate, while for t > 35 min and k < k c it takes a constant 
/c-independent value, interface portions separated by dis- 
tances r > l/k c being uncorrelated. For a rough inter- 
face pL the value of k c decreases with time as k c ~ t~ x l z , 
where z is the dynamic exponent. Indeed, from the log- 
log plot of k c vs t (Fig. 3), a slope value z = 1.6 ± 0.1 
is obtained. Therefore, the time behavior of our CVD 
growth process is well described by that of a self-affinc 
rough surface. The same conclusion is obtained from the 
spatial behavior of the surface morphology. Specifically, 
for a two-dimensional rough interface the PSD behaves 
as Q PSD(fc) ~ l/k 2a+2 , with a the roughness exponent. 
In Fig. 2, as time proceeds, up to three spatial regions 
can be observed: 

i) For t < 50 min we can distinguish two regions: Re- 
gion I, for k > ko, features aj = 0.99 ± 0.04; ko ini- 
tially decreases from 0.04 nm -1 for t = 20 min to a con- 
stant value of 0.02 nm -1 for t > 50 min. Region II with 
an = 0.75 ± 0.03 is observed for k\ < k < ko, where k\ 
changes from 0.025 nm -1 for t = 20 min to 0.0034 nm -1 
for t = 50 min. 

ii) For t > 50 min a new region III appears for 
k c < k < ki, with am — 0.42 ± 0.03; as noted above, k c 
decreases with time. Regions I and II are still observed. 

The above behavior of the PSD is compatible with 
that of the surface width, W 2 {t) = {{h(r,t) -h(t)) 2 ), 
which we measure independently (Fig. 3). The bar de- 
notes spatial average. For a rough interface, initially 
W(t) ~ t a ' z — t 13 , as long as t <C L z , with L a measure 
of the spatial extent of the system. For t ^> L z , the inter- 
face saturates into a stationary state for which W ~ L a . 
As a difference with the PSD, we observe only two tem- 
poral behaviors in W(t), corresponding to regions II and 
III of the PSD. For t < 50 min, the width is described 
by the growth exponent (in = 0.42 ± 0.04, while for 50 
min < t < 15-30 hours, W(t) data are consistent with 
Pin = ctin/z = 0.26 ± 0.03 as obtained from the PSD 
data. For t > 15-30 hours, the width W(t) saturates to a 
constant value. Note that the contribution of the short 
length scales (region I) to the surface width is masked 
by the spatial region II. However, despite the few data 
points for t < 50 min, an attempt to estimate Pi from 
the PSD curves leads to Pi = 0.28 ±0.09. Regarding 
the long distance properties of our system, the observed 



a m = 0.42 ± 0.03, z = 1.6 ± 0.1, and p m = 0.26 ± 0.03 
are close, within experimental errors, to the approximate 
exponents for the KPZ universality class in three di- 
mensions, namely olkpz = 0.39, zkpz — 1-63, and 
Pkpz = 0.24. Thus, the asymptotic scaling behavior of 
the growing CVD surface is of the KPZ class, an experi- 
mental result very rarely found in three dimensions 0. 

It is important to understand the origin of the KPZ 
regime under our experimental conditions. It is known 
p"l| that, for CVD systems with a low surface reaction ki- 
netics (low sticking probability s), lateral growth, which 
is the fingerprint of KPZ behavior, is promoted leading to 
the so-called conformal growth. The mechanism induc- 
ing conformal growth is the reemission of the precursor 
{i.e., depositing) species due to a low sticking probability. 
For our experimental conditions s has actually been re- 
ported JlJ] to decrease with temperature from s = 0.5 at 
573 K down to s = 0.08 at 723 K. Also, the growth con- 
formality of submicron trenches by Si02 films has been 
reported Jl2|,[l3| to improve as T increases. Physically, 
s can be related to the concentration of reactive sites at 
the growing film surface, which for our system are mainly 
associated with hydrogenated groups, such as Si-OH and 
Si-H, and with strained siloxane (Si-O-Si) groups, rather 
than with the less reactive relaxed siloxane groups jl4| . 

To study the role played by the sticking probability 
on our film morphology, we have grown SiC>2 films at 
lower T = 611 K (higher s) with the same growth rate 
(w 20 nm/min). Analysis of our films by Raman and 
infrared spectroscopies shows jll| that the 611 K films 
indeed present a significantly higher concentration of re- 
active groups than the 723 K films, in agreement with the 
higher sticking probability expected for the low temper- 
ature growth. Moreover, regarding the surface dynamics 
of the low T films, neither the KPZ nor the saturation 
regimes were found after 2 days of deposition. In Fig. 4 
we plot the PSD of a film grown at 611 K after 2 days. 
The KPZ region does not appear but only regions I and 
II do. In Fig. 4 (inset) we plot W(t) for the same set 
of films. Again, only unstable growth is obtained for 
long times. From the point of view of Eq. (Q), these ex- 
periments suggest that when T is reduced the effective 
A coefficient of the KPZ term decreases. These results 
are analogous to the crossover found between diffusion- 
limited aggregation and Eden growth when tuning the 
sticking probability fjq| , and allow us to identify region 
II in the scaling behavior at the high tempererature value 
(T = 723 K) as an unstable transient. 

The above results suggest that the KPZ equation might 
be a good starting point to describe the observed behav- 
ior. However, in our system we expect the first term in 
equation ([!]) to be negligible, since the Si02 vapor pres- 
sure is extremely low (i.e. v is very small) in our tem- 
perature range, whereas surface stabilization by surface 
diffusion seems to be relevant ||. In fact, ai ~ 0.99 
and (3i k, 0.28 are consistent within experimental errors 
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with the linear theory of surface diffusion |L7||. On the 
other hand, due to the large values of the effective expo- 
nents an and /3/j and in view of the discussion above, we 
believe region II corresponds to unstable growth. This 
instability can be related to shadowing effects which oc- 
cur in CVD |^-^| due to the random walk motion of the 
depositing particles. Nonlinear surface diffusion pl| can 
be discarded as the origin of the unstable behavior in re- 
gion II since it leads to anomalous scaling 0] , not present 
in our measurements. Similarly, the unstable equation 
for h derived in |6| can be ruled out as a description of 
the unstable behavior in region II because it leads to a 
non constant growth rate, incompatible with our experi- 
mental setup. Thus, we propose the following continuum 
equation to describe the silica CVD growth: 



d t h 



-KV 4 h + S0/9+ (A/2)(V/i) 2 + r)(r, t). (2) 



In Eq. (g), K and e are positive constants; the first term 
on the rhs represents relaxation by surface diffusion and 
the second term represents geometric shadowing effects, 
9 being the local exposure angle and being the spatial 
average of 9 Jl^|. In order to check the validity of (||) 
we have performed numerical simulations in two dimen- 
sions (three dimensional simulations are limited to small 
system sizes [Q). In Fig. 5 we plot the time evolution 
of the PSD for K = D = 1, e = 1/2, and two values 
of the strength of the KPZ non-linearity A = 0.2,3. For 
these sets of parameters we have employed system sizes 
up to L — 8192 in order to reliably identify the different 
scaling regimes. For A = 3 (solid line in Fig. 5), and as 
experimentally observed at high T (Fig. 2), the surface 
PSD evolves from a scaling regime dominated by sur- 
face diffusion at short distances to KPZ scaling at large 
length scales, through an unstable transient due to the 
shadowing effects pi]] . Also in agreement with the ex- 
perimental data, after a certain growth time both ko and 
k\ are frozen, and no anomalous scaling is observed. The 
time evolution of the surface width is shown in the inset 
of Fig. 5 for the same parameters. Again, for A = 3, it 
resembles quite closely that shown in Fig. 3 |2l[] since, 
as t increases, the behavior changes from surface diffu- 
sion to unstable growth and finally to KPZ scaling. For 
smaller A = 0.2 (dashed lines in Fig. 5), corresponding 
experimentally to low T (Fig. 4), we observe from the 
PSD and W(t) behaviors that the crossover from the un- 
stable region II to KPZ behavior takes place, if at all, 
at longer time and larger length scales. Similar behav- 
iors to those in Fig. 5 are obtained for other parameter 
sets 15 1, provided that the relative weight of the different 
contributions is preserved. 

In summary, we have found that low pressure CVD 
growth of silica films is governed by the relative balance 
between surface diffusion, shadowing and lateral (KPZ) 
growth. Our experimental system is well described by a 
differential continuum stochastic equation in which these 



three mechanisms are allowed for. Moreover, our study 
allows to link the value of the effective KPZ nonlinear- 
ity to the physical and chemical properties of the grow- 
ing interface. Thus, we conclude that the observation of 
asymptotic KPZ scaling is favoured under growth con- 
ditions (e.g., high temperatures) which promote a low 
sticking probability of the depositing species. 
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tions. However, the qualitative behaviors are the same: 
for Fig. 5 we obtain for the A = 3 case af = 1.5 ± 0.1 
(/3f = 0.37 ± 0.01), a a n = 0.86 ± 0.04 (/% = 0.45 ± 0.02), 
and am = 0.5 ± 0.01 {f3j n = 0.33 ± 0.02), consistent, re- 
spectively, with linear surface diffusion, unstable growth, 
and KPZ scaling in two dimensions. 




FIG. 1. AFM 3x3 /im 2 image of amorphous SiC>2 films 
deposited at T = 723 K for t — 75 min. The vertical bar 
represents 400 nm. 
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FIG. 2. PSD curves from AFM images of the Si02 surface 
at different times: substrate (solid line), 5 min (x), 20 min 
(•), 50 min (V), 180 min (♦), 480 min (a), 900 min (+), 1800 
min (*), and 2880 min (□). Wave vectors ko, ki, and k c sep- 
arating the scaling regions are shown for t = 2880 min. The 
slopes of the dashed lines are ctni = 0.42 (for k c < k < ki), 
a n = 0.75 (for fei < k < ko), and ai = 0.99 (for k > ko). 
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FIG. 3. Plot of W{t) (■, left axis) and k c (t) (o, right axis). 
The slopes of the solid and dashed lines are, respectively, 
(3m = 0.26 and 1/z = 0.62. 
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FIG. 4. PSD of a Si0 2 film grown for 2 days at 611 K. 
Slopes of the dashed lines are an = 0.58, ai = 0.9. Inset: 
W(t) for the same system. Slope of the solid line is (3n = 0.51. 
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FIG. 5. PSD from numerical integrations of Eq. (Q) for 
K — D — l,e — 1/2, L = 1024, t = 100, averaged over 250 
realizations, for A = 0.2 (dashed line) and A = 3 (solid line). 
Slope of the straight solid line is a s m = 0.5. For A = 3, the 
values af = 1.5, afj = 0.86 give the best fit to the PSD be- 
havior in each region. Data curvature at large k is due to 
discretization effects. Inset: W(t) for the same parameters. 
Slopes of the straight solid lines are /3| = 0.37, /3|j = 0.45, 
and Pfn = 0.33. 
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